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Abstract. A rigorous mathematical proof is given of a class of vector identities 
that provide a way to separate an arbitrary vector field (over a linear space) into 
the sum of a radial (i.e., pointing toward the radial unit vector) vector field, minus the 
divergence of a tensor plus the curl of an axial vector. Such a separation is applied to 
the representation of electric current densities yielding a specific form of the effective 
polarization and magnetization fields which is also discussed in some details. 



On dealing with the description of the electromagnetic response of material media, 
the issue of the representation of the induced current density appears to be crucial for 
the definition of the effective polarization and magnetization fields and, thus, for writing 
the macroscopic Maxwell's equations [Vinogradov and Aivazyan 1999, Vinogradov 2002, 
Agranovich et al 2004, Bornatici and Maj 2005]. Specifically, one should find the 
effective polarization and magnetization fields. Peg and Meg, as functionals of electric 
the current density J so that the basic relationship [Jackson 1999] 



is satisfied (Gaussian units are used). In particular, on accounting for the charge 
continuity equation, the representation ([T]) yields the corresponding relationship for the 
charge density p(r, t), namely, p(r, t) = —V • Pefr(r, t), the required integration constant 
being properly chosen. 

It is readily verified that the effective fields Pofr and Meg verifying equation ([1]) for 
a given current density J are not uniquely determined and thus different forms of the 
effective polarization and magnetization fields are available, obtained on the basis of 
either the Helmholtz's theorem [Bornatici and Maj 2005] or vector identities proved by 
Vinogradov and Aivazyan [Vinogradov and Aivazyan 1999, Vinogradov 2002] and by 
Raab (quoted in [Vinogradov and Aivazyan 1999]). 

In this letter, a rigorous mathematical proof is given of a class of identities which 
are satisfied by any (sufficiently regular) vector field over a linear space. When applied 
to the current density J(r, t) regarded as a vector field over the three-dimensional r- 
space (with time t treated as a parameter), each one of such identities (labelled by an 
integer n > 1) can be cast in the form and, thus, yields a specific form of effective 




(1) 
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polarization and magnetization fields. The results reported by Vinogradov and Aivazyan 
[Vinogradov and Aivazian 1999, Vinogradov 2002] are found as particular cases for n = 1 
and n = 2. 

With this aim, let us consider a generic vector field R(x) = {Ri{x.))f^^ with 
components over the iV- dimensional (linear) space with x = (xi, . . . ,xn) (it is crucial 
that the dimension of the vector field is the same as that of the considered space) and 
assume that the field is sufficiently regular, i.e., different iable up to a high-enough order. 
We claim that 

d d dRk 

rjt , rj^ , , , , . • • • 

' dxj^ dxj„^, dxk 

where repeated indices are summed, n > 1 is an integer, ("^^) denotes the binomial 
coefficient and the following quantities have been defined as functionals of R(x), viz., 

d 



and 



[7^^"^(x,R(x)) = __Q^^ {x^Xk, ■ ■ ■ Xk„^,RkMj , {n>l), (3a) 
gf (x,R(x))=i?,(x), 

Qi (x,R(x)) = — gifc,...fc„(x,R(x)), {n > 1), 



with 



toi-fc„(x, R(x)) = X(iXk^ ■ ■ ■Xfc„_,i?fc„)(x) 

= (^^1)1 ^^'ri^)^'^(H) ■ ■ ■a;<x{fc„_i)i?a(fc„)(x) 

the symmetric part of the tensor XiX^^ ■ ■ ■ x^^^^Rkni'^), the index a running over the 
permutations of the n + 1 indices {i, ki, . . . , kn}. For the case n = 1, equation ([2]) 
reduces to 

dRk _ djxiRk) _ 

OXk OXk 

which can be readily verified. Hence, let us prove the generic case by induction, that is, 
we need to prove that, if ([2]) is true for n, then it is true for n + 1. With this aim, it is 
convenient to consider the identity 



d d dRk d / d d dR, 



dxj^ dxj„ dxk dxj \ ^ -^"^^ dxj^ dxj^_^ dxk 

which, on inserting equation (j2D, yields 
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If • If • • • • T* ■ 



d 



d dRk 



dxj^ dxj^ dxk 



n-l 



i=0 



n-l\ (A^ + n-1)! d 
i ) {N + n- i-l)\dxj 



n-l 



i=0 



n — 1 



{N + ny. 



i J{N + n-i-l)\ 



\J\ ' -{n- l)q\ 



(4) 



Moreover one has 
d 



dxj 
d 



XiU, 



(n) 



where the identity 



(n)\ _ n + 2 



(n) 



n+ 1 

Qn 



n+1 



(5) 



n+1 dxfc, ■ ■ ■ dxk„ \ J 

has been used. Equation (jlj), along with ([5]), reduces to ([2]) with n ^ n + 1, hence one 
can conclude that the identity ^ holds true for any integer n>l. 

By inspection of ([2]) one should note that it is convenient to separate the symmetric 
part of the tensor XiXk^ ■ ■ ■ Xk„_^Rkn in the quantity U^^\ namely, (cf. equation (l3aj) ) 

d 



/7i")(x,R(x)) = 

Qn 

dxk^ ■ ■ ■ dxk„ 



■ — I XjX/jj^ ■ ■ ■ Xj^ ^^Rj^ -\- Xj^X}^^ • • • Rkn—i'^kn ~l~ ■ ■ ■ ~l~ Ri-^ki ' ' ' ■^kn ) 
n + 1 V / 



n + 1 

(") , 

Ti + 



n 1 

•^i-^ki ' ' ' ■^kn—l-^kn ! T ( •^i'^kl ' ' ' Rkn-l'^kn + ■ ■ ■ + Ri^kx ' ' ' ■^kn 



dXkx ■ ■ ■ dXk„ 



n + 1 



n+1 



^X{Rk^ Xk„Ri)Xki ' ' ' Xk„_i 



(n) , d („) 

j 



(7) 



where the vector ql"'' = g|"'^(x, R(x)) is defined in equations (l3&j) and the axial vector 
m^"^ = m^"^(x, R(x)) is given by 

Qn— 1 

<(x.R(x)) = — 

with 



dxki ■ ■ ■ dxk^_^ 



"^feA:i-fe„_i(x, R(x)), (n > 1), (8a) 



"^fcfci-fc„_i(x,R(x)) = —^(ekabXaRbi?^)]xkx " ' ' Xk^.^- 



lb) 



In equation ([7]), the identities ([6]) and eijkSkab = ^ia^jb — ^ib^ja have been accounted for, 
Eijk and 6ij being, respectively, the completely anti-symmetric permutation symbol and 
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the Kronecker delta. On substituting 
the same form, one has 



into ([2]) and grouping together the terms of 



d 



n-1 



d dRk (n) 



9xj, dxj^^, dxk 

^ {N + n-l)\ fn - 1 
^ (A^ + n-^ - 1)! V^ - 1 



n 



N 



N + n 



_d_ 

dxi 



m 



in) 



m 



in-i) 



+ ^"^^ m 



(9) 



where equations fl3B for n = has been used. Equation ^ can be solved for R,, as a 
function of its own derivatives. To this end it is convenient to change the dummy index 
I into n — I with the result that, for any integer n > 1, 

R = X(")-V-T(") + Vx Y("\ (10) 

where the fields X(") = T^") = (T4"^)ffe=i and Y^") = are functional 

of R(x) given by 



xf")(x,R(x)) 



TV! 



7;i")(x,R(x)) = 



d 



d dR,U) 



■X 



n 



{N + i-l)\ \i 



^" ^dxj, dxj„_, dxk 
{n-e){N-l) 



[lla) 



dxki ■ ■ ■ dxk,_. 



m 



n{N + i) 
gifcfcr..fc,,i(x,R(x)), 



{N + ndxk,---dx 



[Ub) 



"lfcfci...A:,_i(x,R(x)). (lie) 



Equation ( fTOl) constitutes the main result of this paper according to which a generic 
vector field R(x) can be decomposed into a radial vector X'^"\ pointing toward the radial 
unit vector x = x/|x| (with |xp = ^^e divergence of a tensor T*^") and the curl 

of an axial vector Y^'^^ It is worth noting that for every integer n > 1, there is a set of 
fields X*^"\ T*^"^ and Y*^"^ which satisfy identically equation (fTUj) . the right-hand-side of 
( ITOl) being, indeed, independent on n. 

The identity (fTOj) can be readily applied to represent the current density J(r, t) over 
the three-dimensional r-space, the time t being treated as a parameter. Specifically, for 
fixed time t, one can set J(r,t) = R(x), with r = x = (xi,a;2,a;3) and = 3, and, for 
any integer n > 1, define the effective polarization, Pg^'', and magnetization, M^^'', by 

»(")/ 



X(")(r, J(r,t)) - V ■ T(")(r, J(r,t)) = 



and 



Y(")(r,J(r,t)) = cMiJ(r,J(r,t)), 

respectively, so that equation (fTUj) amounts to the representation ([T]) for the electric 
current density J(r, t). 
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On accounting for equations (IllfflP -i nTcj) explicitly, the effective fields P^^'' and 

In) 

can be expressed in a form which, formally, resembles the multipole expansion 
[Agranovich et al 2004], namely, 

n-ff .(r,t) = P/"Hr, J(r,t)) -— E ^.. Qtl£-..-S-^^i-^^^ (12a) 



dxk <9xfci ■ ■ ■ dxki_^ 



[I2h) 



where Qikki---ki_i ^kk'i--ke_iP^^y '^o^^ of electric dipole, electric quadrupole 
and magnetic multipole moments, respectively, and are functionals of the induced 
current density J(r, t) given by 

dP'r^, / 6 d d dJkirA) 



dt 

An/) 



[r,J(r,t)) = (-l) 



6 In 



2(n-f) 
n(£ + 3) 



(£ + 2)! \i 
6 /n\ £ , J ( w 
c(£ + 2)! y'n(£ + l)^^'^'^^ ^^^^'^ ' ' 



a;(ia;fci ■ ■ ■ x^^.i Jfc)(r, t), (136) 

(13 c) 



Equation ([T]) along with the effective fields (12) takes the form 



_dP^_d__d_^ 



^ ikki---ki_i 



dt dt dxk Sxfci ■ ■ ■ 



dxi 



d 



,1-1 



'dxj 
J 1=1 



dxk.---dxk, , 



(14) 



which exhibits the same mathematical structure as a multipole expansion truncated to 
the n-th order, even though, with the "multipoles" given by equations fll3fflp - fll3cp . it is 
an exact identity. 

For the specific case n = 1, equation (fT^ reduces to 



J,; 



dJk d 



~Xi 



dxh dxi 



d 



dx 



J L 



2c 



^ kmnXmJn 



(15) 



which is the identity proved by Vinogradov and Aivazyan [Vinogradov and Aivazyan 
1999]. 

For n = 2, equation ffT^ reads 



(2,1) 



(2,1) 



dt dt dxj dt dxjdxk 



dXj 



^2^(2,2) 



with 



dP 



(2) 



1 



d dJk 



dt 



(r,J(r,t)) = -x.x,-^ 



dQp^^ 3 

(r,J(r,t)) = --{xiJj + XjJi 



(17a) 
(176) 
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(r,J(r,t)) = — (xjXjJfc + XiJjXfc + JiXjXk), (17c) 

iv4'''^(r, J(r,t)) = —{eij^XjJk)xi. {lie) 

Equation (fT6l) . along with (17), agrees with Raab's result quoted in [Vinogradov and 
Aivazyan 1999]. (The factor 1/2 in the expression ffTTdl) for M(2'i) appears to be missing 
in the corresponding Raab's expression.) 

With reference to fll3ap and (1136|) , only the time-derivative of both the vector P-''"^ 
and the tensors Q^iktl--ki_i specified as functionals of the induced current density 
J(r,t). The explicit form of both P/"^ and Qikki---ki_i obtained on integrating with 
respect to time, 

dt'^{r,J{r,t')) + Pt\r), (18) 

QLt....._,(r, J(r,t)) = dt' ^^^^;;-^--H r,J(r,tO) + gLt...,,_,(r), (19) 

where to is a reference time (as it appears below, its presence is crucial for the treatment 
of static fields) and /"/"^r) and Qikki---ke_i{^) time-independent integration constants 
which should be selected in such a way that the representation in ([T]) for the induced 
charge density p(r, t) is verified, i.e., V ■ P^^{r,t) = —p{r,t). 

In particular, the integration in ( |T8l) can be explicitly carried out on making use of 
fll3ap along with the charge continuity equation, with the result that Pj^"^ is a functional 
of the induced charge density, 

6 d d 

+ ^/"V). (20) 
For the specific case of fluctuating perturbations, it is usually assumed that both p(r, t) 
and J(r, t) tend to zero for t -co, so that one can set to = -co in equations (fTSj) and 
( |T9|) as well as ( l20l) . with both P/"''(r) and Qikki---ki_ii^) being zero. On the other hand, 
for the case of static fields, for which, in particular, p(r,t) = p(r,to) = p(r), equation 
(1201) reduces to 

p/")(r,p(r)) = p("Hr), (21) 
and equation (fT9ll along with l\13bh where J(r, t) = J(r) yields 

+ QtkiUj^)- (22) 
One should note that equation ( l22l) for the effective electric multipoles exhibits a time- 
dependence even though static fields are considered. As for the effective magnetization, 
the static case is simply given by fllScp with J(r, t) replaced by J(r). 



X{iXk^ ■ ■ ■Xk,_^Jk){r) 
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Let us examine in some detail the effective fields fll2ap and (11261) for the specific 
case of n = 1 [Vinogradov and Aivazyan 1999]. With n = 1, the effective polarization 
( IM is 

P^J^{r,t) = P«(r,p(r,t)) - V ■ Q^'''\r, 3{r,t)), 
with, from (I^Uj). 

p(r,t)r, 



(23) 



P«(r,p(r,t)): 
and, from (fT9l) along with ( 113&I) , 



(24) 



p(i)(i 



I' rft'(rJ(r,t')+J(r,t') 



Q(i'i)(r,J(r,t))= <^ 



(25) 



(t - to) 



(rJ(r) + J(r)r) +Q(i'i)( 



the upper (lower) entry of fl^^ and fl2^ referring to fluctuating (static) fields. On 
making use of the charge continuity equation, from ( l25i) one gets 



V-Q(^'^)(r,J(r,t)) = < 



ip(r,t)r- 1* rft'(2J(r,tO + ^(r-V)J(r,t')): 
- '^o) ^^j^^^ ^ _ V)J(r)) + V ■ Q(^'^)(r). 



(26) 



It is readily verified that the effective polarization fl23p along with ([21]) and 0251) is such 
that V-Pj,J:V>^) = -P(r,t) this result requiring that V- (P(^)(r) - V- Q(^'^)(r)) = -p(r) 
for the case of static fields. 

As for the magnetization ( 1126|) . for n = 1, it is 



1 



MS(r,t) = M(i'i)(r, J(r,t)) = -(r x J(r,t)), 



(27) 



which agrees with the conventional definition of magnetic moment density associated to 
the current density J(r, t) [Jackson 1999, Vinogradov and Aivazyan 1999]. 



References 

Bornatici M and Maj O 2006 Physica Scripta 73, 160 
Jackson J D 1999 Classical Electrodynamics (New York: Wiley) 
Vinogradov A P and Aivazyan A V 1999 Phys. Rev. E 60, 987 
Vinogradov A P 2002 Physics Uspekhi 45, 331 



